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Clock synchronization between the ground and satellites is a fundamental issue in future quantum telecom¬ 
munication, navigation, and global positioning systems. Here, we propose a scheme of near-Earth orbit satellite- 
based quantum clock synchronization with atmospheric dispersion cancellation by taking into account the space- 
time background of the Earth. Two frequency entangled pulses are employed to synchronize two clocks, one 
at a ground station and the other at a satellite. The time discrepancy of the two clocks is introduced into the 
pulses by moving mirrors and is extracted by measuring the coincidence rate of the pulses in the interferometer. 
We find that the pulses are distorted due to effects of gravity when they propagate between the Earth and the 
satellite, resulting in remarkably affected coincidence rates. We also find that the precision of the clock synchro¬ 
nization is sensitive to the source parameters and the altitude of the satellite. The scheme provides a solution 
for satellite-based quantum clock synchronization with high precision, which can be realized, in principle, with 
current technology. 
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I. INTRODUCTION 

High precision synchronization of clocks plays an impor¬ 
tant role in modern society and scientific research Cl 11 ; ex¬ 
amples include navigation, global positioning, tests of general 
relativity theory, long baseline interferometry in radio astron¬ 
omy, as well as gravitational wave observation. Two standard 
classical protocols for clock synchronization are Einstein’s 
synchronization scheme |H and Eddington’s slow clock trans¬ 
fer @]. The former requires operational exchange of light 
pulses between the distant clocks and the latter is based on 
sending a locally synchronized clock from one part to other 
parts. Recently, quantum strategies have been exploited to 
improve the accuracy of clock synchronization. A few quan¬ 
tum clock synchronization (QCS) proposals JUJU and exper¬ 
iments are reported. It is shown that the schemes based 
on quantum mechanics can gain significant improvements in 
precision over their classical counterparts. 

On one hand several satellite-based quantum optics experi¬ 
ences M are feasible with current technology, such as satel¬ 
lite quantum communication fl9l - l27ll . and quantum tagging 
[28], as well as gravity probes using beam interferometers 
129, |30] and atomic clocks sum to test the principle of 
equivalence. Among these experiments, the synchronization 
of clocks between a satellite and a ground station |33| is an 
essential step. Besides, a satellite-based quantum network of 
clocks is promising to act as a single world clock with un¬ 
precedented stability and accuracy approaching the limit set 
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by quantum mechanics, and there is also a security advan- 
tage 0. On the other hand, time dilation |32 j| is a concern 
because of relativistic effects of the Earth on the QCS. The 
influence of relativistic effects on quantum systems 05l - t43[l is 
a focus of study in recent years because such studies provide 
insights into some key questions in quantum mechanics and 
relativity, such as nonlocality, causality, and the information 
paradox of black holes. Relativistic effects are particularly 
significant for the quantum versions of the Eddington scheme 
[6H8[] because one must assume that the transfer is performed 
“adiabatically slowly” [44]] and the spacetime is flat such that 
relativistic effects are negligible. However, time dilation in¬ 
duced by Earth’s spacetime curvature is experimentally ob¬ 
served for a change in height of 0.33m [32] and thus cannot 
be neglected. 

In this paper we propose a practical scheme for satellite- 
based QCS. We let two observers, Alice and Bob, exchange 
two frequency entangled pulses between a ground station and 
a satellite. The influence of gravitational red-shift on the fre¬ 
quency of a pulse can be eliminated by an opposite gravi¬ 
tational blue-shift. By assuming the clocks have the same 
precision EM, clock synchronization can be realized by 
identifying the time discrepancies. In our scheme, the time 
discrepancies are introduced through adding or subtracting 
optical path differences into the pulses and the coincidence 
rate of the pulses in the interferometer as a function of the 
time discrepancy between the clocks is measured. In actual 
satellite-based quantum information processing tasks El [Hi- 
Ill 111 min, and similarly in protocols of QCS, the main 
errors are induced by photon-loss and the dispersion effects 
of the atmosphere through which the pulses travel. Therefore, 
we employ frequency entangled light, instead of the entangled 
NOON state which is vulnerable for photon-loss J0, as well 
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as the dispersion cancellation technology fsl 1471 l48ll to elimi¬ 
nate the influence of the atmospheric scattering. We find that 
the coincidence rate of interferometry is remarkably affected 
by the spacetime curvature of the Earth. We also find that the 
precision of the clock synchronization is sensitive to the light 
source parameters. 

The outline of the paper is as follows. In Sec. II we briefly 
introduce the sketch of the experimental setup. In Sec. Ill we 
discuss how the Earth’s spacetime effects the propagation of 
photons. In Sec. IV we study the feasibility of satellite-based 
QCS and how the effect of the Earth’s gravity will disturb it. 
In the last section we discuss the experimental feasibility of 
our scheme and give a brief summary. 

II. THE SCHEME 

The sketch of our proposal for satellite-based QCS is de¬ 
scribed in Fig. (|T|). The quantum optical technology of our 
proposal is based on the Hong, Ou, and Mandel (HOM) in¬ 
terferometer jH. We assume that Alice works on the sur¬ 
face of the Earth (r = r a) with her own clock, while Bob 
works on a satellite at constant radius r = vb > r a- The 
clocks have the same accuracy thus the clock synchroniza¬ 
tion problem is reduced to the problem of identifying the time 
discrepancy between the clocks. Alice sends two frequency 
entangled beams produced by a parametric down converter 
crystal (PDC) source to Bob and Bob bounces them back to 
Alice again. Those two pulses are named signal ( S ) beam 
and idler (/) beam, respectively. By exchangingentangled 
beams between Alice and Bob, a “conveyor belt” 131 for time 
information is established. After propagating through differ¬ 
ent optical paths, the signal and idler beams are interfered at 
the 50/50 beam splitter and be measured by the detectors. To 
introduce time information into the beams, Alice and Bob use 
moving mirrors with constant speed v to add or subtract op¬ 
tical path differences (OPD) to the beams @]. Alice and Bob 
come to an agreement on the starting time To of their mirrors 
in advance. Since they do not have a synchronized clock to 
start with, they can only start the mirrors at time Tq(* = a, b) 
relative to the time readings of their own clock, which are dif¬ 
ferent due to different locations. As described in Fig. (1), Al¬ 
ice at the ground station start a moving mirror to add an OPD 
Sl^(t) to the idler beam fU and to subtract the same amount 
of OPD 51% (f) from the signal beam which was bounced back 
from Bob. At the satellite Bob subtracts an OPD 6ll(t) from 
the idler beam and adds an identical OPD 51% (t) to the signal 
beam. The linear time dependent OPDs are given by 

K(t) = v(t - Tg ), Sll(t) = -v(t - Tq ), 

( 1 ) 

M%(t) = ~v(t - Tq ) , 61$ (t) = v(t - Tq ) , 

where Tq and Tq are the starting times (proper time) as mea¬ 
sured by Alice’s and Bob’s clocks, respectively. In Eq. 0} t is 
the time point (coordinate time) of the coincidence detection 
of the signal and idler photons, i.e., the time when the pho¬ 
tons quantum state is measured. By assuming the quantum 
state is instantaneously collapsed by the measurements made 


on the surface of the Earth, we can agree that the collapsing 
time is identical t = t a = t b , even if Alice and Bob’ reading 
times are different r“ / r b at this moment. From Eq. (|7]i 
we can see that the delays are proportional to the time inter¬ 
val between the starting time Tq for the moving of the mirrors 
and the time t when the photons are detected. If the propor¬ 
tionality constant v in Eq. Q and the starting time reading 
Tq on Alice’s and Bob’s clocks are identical, the quantity of 
OPD at Alice’s point will be zero after an exchange period. 
Therefore, the final OPD is totally produced from the starting 
time discrepancy At = Tq — Tq between Alice’s clock and 
Bob’s clock. Then the signal and idler pulses are interfered at 
the beam splitter (BS) and click the detectors. We will show 
in section IV that the final difference of optical path lengths 
is affected by a factor depending on the starting time discrep¬ 
ancy At. By measuring the photon coincidence rate P c at the 
output ports 1 and 2 of the beam splitter, one may acquire very 
precise information on the OPD in the two arms. Thus, it is 
sufficient to measure the photon coincidence rate to recover 
the exact time discrepancy between Alice’s clock and Bob’s 
clock. Then Alice tells Bob by classical communication to 
adjust his clock according to the time discrepancy. By using 
this scheme, Alice and Bob may check how much the accu¬ 
racy is disturbed by the gravity induced spacetime curvature 
of the Earth, and inversely they can precisely measure the cur¬ 
vature via an atomic clock, which is the most accurate setup 
currently available in the world. 

HI. EARTH’S SPACETIME ON THE PROPAGATION OF 
PHOTONS 

Now we describe the propagation of photons from the 
Earth to a satellite by taking gravity of the Earth into con¬ 
sideration. The Earth’s spacetime curvature will influence 
the light pulses during their propagation between the ground 
station and the satellite. We know that the Earth rotates 
slowly with an angular velocity at the equator of oje = 
7.2921151247 x 10 -5 rad/s or at a linear speed of ve = 
465 ms -1 , which is much slower than the speed of light. 
Therefore, the Schwarzschild metric toj - lvil l is a sufficient ap¬ 
proximation for the Earth’s spacetime, as has been discussed 
in lf27l [45l [49il . The Schwarzschild metric is given by lf50ll5lll 

d S 2 = _( 1 _ 11)^2 + ( 1 _! l ) - ldr 2 
r r 

+r 2 (dd 2 + sin 2 0dtp 2 ), (2) 

where rg := 2< ^ M is the Earth’s Schwarzschild radius, M is 
the mass of the Earth, c is the speed of light in vacuo, and G 
is the gravitational constant. 

A photon can be properly modeled by a wave packet of 
electromagnetic fields ll53h with a distribution F^) 0 of modes 
peaked around the frequencies uik ,0 f54L 56], where K = 
A, B labels either Alice or Bob. The annihilation operator 
for a photon from the point of view of Alice or Bob takes the 
form 

r+00 

<W. 0 (t K )= ckvft K e- iUKtK F^ o (uj K )a UK , (3) 

Jo 
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FIG. 1: Sketch of the experimental setup for satellite-based QCS. Alice sends two frequency entangled beams produced by a parametric down 
converter crystal (PDC) source to Bob who bounces them back to Alice again. Alice adds an OPD Sl^(t) to the idler beam and subtracts 
the same amount of OPD (t) from the signal beam which was bounced back from Bob. At the satellite Bob subtracts an OPD <5 1{ (f) from 
the idler beam and adds an identical OPD (f) to the signal beam. Finally, we let the pulses interfere at the 50/50 beam splitter (BS) and be 
detected at the measurement setups. The coincidence rate P c is measured as a function of the time discrepancy tq — t§ between the clocks. 


where lok is the physical frequencies as measured in their 
labs. At time ta Alice prepares a wave packet F^ A ' Q and 
sends it to Bob who will receive it at a different proper time 
t b = Ar + \Jf (r B )/f (r a )ta, where f(r A ) = 1 - 
At is propagation time of the wave packet, and the factor 
yj firs')/f (t'a)'ta indicates the relativistic time dilation. The 
wave packet received is modified due to the spacetime cur¬ 
vature of the Earth and takes the form F ^ 0 . By employing 
the definition of proper time, it is easy to show that the time 
evolution for the modes has the form id TK (j) u]K = 

(u) 

where <f>u£ are the quantum states of the modes correspond¬ 
ing to the operators a UK 0 This equation indicates that the 
physical frequency u>k measured by an observer at the posi¬ 
tion tk is uik = f(fK ) _1 / 2 w. Then we find that Bob will 

receive a mode with frequency ojb = \It$A w a if a sharp 
frequency mode with u>a was sent by Alice. This is the well- 
known gravitational red-shift effect, which was predicted by 
Einstein in 1911 and experimentally verified in 1960 [55}]. In 
our scheme such a classical gravitational effect is designed 
to be eliminated by an opposite gravitational blueshift fac¬ 
tor yj because the signal will be sent downward from 
a satellite to the Earth. 

However, such a nonlinear gravitational effect is found to 
influence the fidelity of the quantum channel between Alice 
and Bob lf27jl45tl and will inevitably affect the accuracy of the 
satellite-based QCS as well. The mode a! received by Bob 
can be decomposed in terms of the mode a prepared by Alice 
and an orthogonal mode a± 10. l45l[57h 

o' = 0a + \/l — 0 2 aj_, (4) 

where 0 is the wave packet overlap between the distributions 

F^b\ ( w b) and F^ 0 (w B ), 

/•Too 

0=/ d^AF^MF^B), (5) 

To 

which describes the fidelity of the channel between Alice and 
Bob. For a perfect channel one has 0 = 1. 


IV. SPACETIME CURVATURE ON SATELLITE-BASED 
QCS 

The emitted signal and idler beams from the PDC initially 
share an entangled state |[s|] 

\ip) = J dcuidLO2F(uJi,uj2)a\u}i)a\uj2)\0), ( 6 ) 

where a'(u> i) and a^u^) are creation operators of the first 
and the second photons, respectively. For a monochromatic 
pump this state can be rewritten as 

| tp) = J du F U0 (u) |w 0 +w)/|w 0 - u) s , (7) 

where \uj)s and \oj)j are the states for the signal and idler 
pulses, and F Uo (w) is the spectral distribution function of the 
down-converted photons |47l l48f |. Now let us briefly dis¬ 
cuss how quantum entanglement is useful to quantum clock 
synchronization and the advantage of the quantum clock syn¬ 
chronization scheme. The main advantage of a quantum strat¬ 
egy is that we can employ quantum uncertainty and coher¬ 
ence time of the frequency entangled photons. From Eq. (7) 
we can see that although the sum frequency 2wo is certain, 
the down-shifted frequencies |cuo — w)s and |u;o + ut)i are 
highly uncertain. The frequencies are largely determined by 
the pass bands of the interference filters that inserted in the 
down-shifted beams ll~52ll . These pass bands have been found 
on the order of 5 x 10 12 Hz, which corresponds to a coherence 
time for each photon on the order of 100 fs [52]. Therefore, 
it is able to measure a time interval of or better than the co¬ 
herence time of the photons (50 fs), with an accuracy of 1 fs 
(10 -15 s). Also, the most important task for the quantum clock 
synchronization is to determine the time interval by measuring 
the rate at which photons are detected in coincidence, which 
relates to the coherence length of the photon wave packet and 
entanglement of the photons. 

As stated before, Alice adds an OPD 6l^(t) to the idler 
beam iU and subtracts the same amount of OPD 61^ (t) from 
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the signal beam which was bounced back from Bob. At the 
satellite Bobsubtracts an OPD SI J b (t ) from the idler beam and 
adds an identical OPD Sl^(t) to the signal beam. Then we let 
the pulses interfere at the 50/50 BS and be measured in the 
HOM interferometer. The detected coincidence rate P c at the 
detectors is given by the Mandel formula |58] 

P c oc J^dhdt 2 (V’|^“ ) ^" ) 4 +) ^i (+ V) , (8) 

where t is the interaction time of the detectors, and the elec¬ 
tromagnetic fields at time tj at the output of the beam splitter 
can be defined by 


Bob and the BS, and 82 is the mode overlap between the dis¬ 
tributions when the pulse is received from the satellite. In Eq. 
([TQi the terms nf and Kj describe the effect of the disper¬ 
sive atmosphere on the signal beam on their way to and from 
the satellite, respectively. Notice that the Doppler shift intro¬ 
duced by the first mirror is to jx and the frequency itself is 
revaluated at to again because the second mirror moves in the 
opposite direction J 8 j]. The analogous procedures can be ap¬ 
plied to the idler beam propagation process, yielding the final 
operator transformations 

a"( w ) = 0 i 0 2 a/(w)e- iajT3+iK ‘ ( “ /x)+iK / (aj) + Aai{u) , 


r< + > = 

= (£H)\forj = l,2, (9) 

where £ is the beam cross section, and x 3 is the position of 
the moving mirrors detector. The relationship between the 
output and input fields of the beam of the moving mirrors can 
be obtained by performing Lorentz transformations on the in¬ 
put fields @j. For example, the factor 
after the beams outflow from the moving mirrors evolves to 
,Jxa 3 ( X to)e~' lulT ^ ( ' T3o ~ x:io ^ c \ where /3 = ^ is the Lorentz 
transformation factor, x = 7^-7 denotes the Doppler shift in¬ 
troduced by the moving mirrors, and x 3 0 is the distance of the 
beam’s delay. For the time delays defined in Eq. (1), we have 
x J 0 = Xu for all j. Then the beams are sent from the Earth 
to the satellite and are influenced by the dissipation of the at¬ 
mosphere and the spacetime curvature. The former produces 
a phase discrepancy in 3 t (to ) and the latter can be described 
by Eq. (5). Taking the signal beam as an example, the full 
procedure can be expressed by the transformations between 
the annihilation operators when the photons are sent from the 
Earth to the satellite 


a's(uj) = &!y/x a s ( X 


/( r ^) A 


f( r B ) 




( 10 ) 


where T 1 = {tq — x$/c) — L/c and 01 is the wave packet 

overlap between the distributions when the pulse is sent from 
the Earth to the satellite. We can see that the annihilation oper¬ 
ator as (to) at Alice’s laboratory evolves into a' s (ui) when ob¬ 
serving at Bob’s laboratory, where L is the distance between 
the Earth and the satellite. Again, the annihilation operator for 
the signal beam evolves into a'g(to) before entering the 50/50 
BS 


«sM 


0i0 2 a/(w)e^ T2+iK ‘ S (“)+«?("/*) 

+ Aa sM . (11) 


wher e T 2 = ^(r 0 ° - r 0 b + ^) + (L/ X + L')/c, A = 
0 2 i/l — 0 ^ + \J 1 — © 2 , L' denotes the distance between 


where T 3 = (tq — Tg — ^ ) — (L jx + L') /c. Because the 

distance between the BS and the PDC source is much smaller 
than the distance between the satellite and the Earth, we as¬ 
sume that L = U . At the output of the 50/50 BS, the modes 
are found to be ai(wi) = -^[ia'j(ooi) e _ *“ l5i / c + ag(wi)] 

and a 2 (to 2 ) = ^[/a's(w 2 ) + a"(u; 2 ) e _l “ 2 5 P C ], where SI is 
the delay introduced to relate the coincidence rate P c with the 
path length. Then the coincidence rate P c defined in © is 
obtained as 


oc / dwic?u; 2 (^|a{(u;i)a 2 (w 2 )ai(u;i)a 2 (u; 2 )|^) 


= / dwidw 2 |(0|ai(wi)a 2 (w 2 )|V’)| 2 . 


( 12 ) 


The matrix element (0|ai(wi)a 2 (o; 2 )|V’) in (IT2l) is given by 


(0|ai(wi)a 2 (u; 2 )|^) = - (0i0 2 ) 2 5(wi + to 2 - 2to 0 ) 


e l{p (j)(toi - too) 


1 — e 


2i(oji —(jJo)'T 4—iAk(uji) 


(13) 


where T 4 = i+^( T o — T o ) — Sl/c, and tp is an overall phase 
term that will disappear by taking the modulus @ 1 , and the 
contribution of the dispersion terms is 


A n(to) = K.f (to) — Kf(to) + Kf(to') — (to') 



where to' = 2to 0 — to. If the properties of the beams propagat¬ 
ing through different optical paths are such that nf = k.j and 
Kj = k{, the dispersion effect of the atmosphere is erased. 
Such conditions can be satisfied by allowing the “from” idler 
beam to propagate at a distance less than the spatial inhomo¬ 
geneities of the atmosphere from the “to” signal beam and, 
equivalently, by allowing the “to” idler beam to propagate less 
than the “from” signal beam @], Here, the “from” beam de¬ 
notes the beam from Alice to Bob and vice versa. Then the 
dispersion suffered by one of the photons can cancel that suf¬ 
fered by the other photon. These two photons can remain 
totally coincident after propagating through different optical 
paths. Substituting Eq. (IT3l > into Eq. (IT2l >. we can obtain 

Pc = (Si 0 2 ) 2 [ dio\P\ 2 (l - cos[—- 61)]), (15) 
J C 1 + /J 
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FIG. 2: (Color online) The relative disturbance of coincidence rate 
A p as a function of the altitude of the satellite tb ■ The parameters of 
the PDC light source are fixed as u>o = 700THz and a = 100MHz. 


FIG. 3: (Color online) The relative disturbance of coincidence rate 
A p as a function of the source parameters u>o and a for the fixed 
distance rg = 6.771 x 10 6 m (LEO). 


where T = F^ 0 (ui). From Eq. (IT5l > we can see that the co¬ 
incidence rate P c is directly related to the time discrepancy 
At = Tq — Tq of Alice’s and Bob’s clocks. That is to say, 
Alice’s and Bob’s clocks can be synchronized by using the 
measured coincidence rate of interferometry. 

To be explicit we next only consider the case in 
which F Uo (u>) is a Gaussian wave packet F^ 0 (w) = 

I _ (^-^o) 2 

I, „ e , where a is the Gaussian width. The wave 

V 27TCT 2 

packet overlaps 0i and 02 are found to be 


©1(2) - 


2A 


1 ( 2 ) 


1 +A? 


'4- 2 [ 1 +^ (2)] 


1 ( 2 ) 


(16) 


where A^) = 1 ± i? and the signs ± occur for ca < rg 

or ta > tb- In Eq- <fl~6t . we define •& = — 1 

and uja,o = w s,o = w o I s assumed. The modes will be 
perfectly overlapped (0 = 1) when Alice and Bob are in 
a flat spacetime /(ra) = /( r s) = 1- For the typical 
resources used in quantum optics experiments, the relation 
$ <C (^ a ) 2 -C 1 should be satisfied ll27l PfsL [59tl . which 

o2 2 

yields 0i = 02 ~ 1-Then we find that the coinci¬ 

dence rate P c has the form 


Pc = (1 - 


^oV ri _ p -2<r 2 (6l-^) 2 /c 2 
8 a 2 ’ 1 


(17) 


We can see that the coincidence rate P c has the factor (1 — 

l9 2 C<7 2 

‘ 8 ff ) 4 compared to that of the flat spacetime case, where 

Pf = 1 - e” 2 ^ {5l ~ w) 2 /° 2 0. We define the effect of 
spacetime curvature on the accuracy of clock synchronization 
as the relative disturbance of coincidence rate 


A p = 


n - p c 

pi 


i- a 


» 9 2 /., 2 
v ^0 \4 

8cr 2 ' 


( 18 ) 


It is now clear that the relative disturbance of coincidence rate 
A p depends on the spacetime parameter •& and the character¬ 
istics of the PDC source. 


In Fig. Q, we plot the relative disturbance of the Earth’s 
spacetime curvature on the coincidence rate A p as a function 
of the distance r is between the satellite and the Earth’s core 
for the fixed light source parameters. It is shown that A p in¬ 
creases as the distance increases, i.e., the accuracy of clock 
synchronization depends on the altitude of the satellite, which 
also verifies that the spacetime curvature remarkably influ¬ 
ences the running of the atomic clocks |[32]. This result is 
very different from that of Ref. 0, in which the coincidence 
rate is independent of the distance L between Alice and Bob 
when the spacetime curvature of the Earth is not considered. 

In Fig. (0, we plot the relative disturbance A p over the 
peak frequency ojq and bandwidth cr of the PDC. It is shown 
that the disturbance of accuracy depends sensitively on the 
bandwidth a of the source, which is similar to the flat space- 
time case 0]. However, here we find that the disturbance on 
accuracy also depends on the peak frequency of the pulses, 
which is different from that of 0 where the accuracy is inde¬ 
pendent of the peak frequency. In this paper we are particu¬ 
larly interested in two typical cases, in which the QCS are per¬ 
formed between the ground station and either a low earth orbit 
satellite (LEO), or a geostationary earth orbit (GEO) satellite, 
respectively. 

The LEO case: The typical distance from the Earth to a 
LEO satellite is about 400km, which yields r,\ = 6.371 x 
10 6 m and re = 6.771 x 10 6 m. Considering that the 
Schwarzschild radius of the Earth is rs = 9mm, it is found 
that <5 ~ “ 7(777 — 77 ) = 4-17 x 10 -11 . We employ a 
typical PDC source with a wavelength of 369.5 nm (corre¬ 
sponding to to 0 = 812THz) and <7 = 100MHz (the relation 
S -C (%^) 2 -C 1 is satisfied). A light source with such 
peak frequency and bandwidth is available, for example, in 
trapped ion experiments [50. The relative disturbance of the 
spacetime curvature on the coincidence rate is obtained as 
A^ = 5.73993 x 10 -8 . The achievement of an optical lattice 
clock with accuracy at the 10” 18 level has been reported in 
Refs. |[60L [gTTi . If we would like to synchronize two clocks up 
to a time discrepancy of r = Ins (as in 0 , or a much lower 
level t = 100 ns as presented in 0 ), the correction of the 
Earth’s spacetime curvature effect will reach 10” 17 s during a 
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single synchronization process. Such a correction is compara¬ 
ble to the accuracy of the atom clocks and thus should be con¬ 
sidered for the QCS between clocks in future satellite-based 
applications. Therefore, we can safely arrive at the conclu¬ 
sion that the spacetime curvature is not negligible when the 
synchronization is performed by LEO satellites. 

The GEO case: The typical distance L between a GEO 
satellite and the ground is about 3.6 x 10 7 m. Therefore, the 
distance between the Earth and the satellite vb is about r/j = 
42.371 x 10 6 m, which yields S ~ — i(— — —) = 6 x 10 -10 . 
In this case the relative disturbance of the spacetime curvature 
on the coincidence rate is A® = 1.18729 x 10 -5 . We find 
that the disturbance of the spacetime curvature on the coinci¬ 
dence rate for the GEO satellites becomes even more remark¬ 
able than that of the LEO case. We remark that the current 
GPS satellites have rs ~ 2.7 x 10 7 m. In the QCS scheme, 
the spacetime curvature is also remarkable. 

Error analysis: It is to mentioning that the velocity vari¬ 
ations of the moving mirrors may induce some errors on the 
coincidence rate. However, the order of magnitudes of the 
movement speed of the mirrors is much smaller than the ve¬ 
locity of light, let alone the velocity variation of the mirrors. 
To be specific, the typical velocity of the moving mirrors is 
10” 1 2 m/s. Let us suppose that this velocity has one percent of 
variation, say 10 -3 m/s, which is much smaller than the veloc¬ 
ity of light. Note that the relative disturbance of the spacetime 
curvature on the coincidence rate is on the order of 10 3 4 5 6 7 8 9 10 for 
the LEO satellites and of the order of 10 -5 for the GEO satel¬ 
lites, which are at least 3 orders of magnitudes larger than that 
of the velocity variations of the mirrors. Therefore, the sys¬ 
tematic errors induced by the velocity variation of the mirrors 
can be safely ignored in the scheme. 

V. DISCUSSIONS 

We have proposed a practical satellite-based QCS scheme 
with the advantages of dispersion cancellation and the ro¬ 


bust frequency entangled pulses of light by taking the effects 
of the spacetime curvature of the Earth into consideration. 
The spacetime background of the Earth is described by the 
Schwarzschild metric, and the quantum optics part of our pro¬ 
posal is based on the HOM interferometer. By eliminating 
the gravitational redshift and blueshift of the laser pulses and 
the atmospheric dispersion cancellation, the accuracy of clock 
synchronization in our quantum scheme can be very high by 
showing that A p is close to unity. Our proposal can be im¬ 
plemented, in principle, with current available technologies. 
To be specific, optical sources with the required peak fre¬ 
quency and bandwidth have been achieved by the trapped ion 
experiments iBoll . The feasibility of photon exchanges be¬ 
tween a satellite and a ground station has been experimen¬ 
tally demonstrated f2lj] by the Matera Laser Ranging Obser¬ 
vatory (MLRO) in Italy. Most recently, they have reported the 
operation of experimental satellite quantum communications 
H22I1 by sending selected satellites laser pulses. Our scheme 
can also be generalized to the quantum clock network cases 
(SHI. The results should be significant both for determining 
the accuracy of clock synchronization and for our general un¬ 
derstanding of time discrepancy in future satellite-based quan¬ 
tum systems. 
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